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Cold atom traps and certain neutron star layers may contain fermions with separation much larger than the 
range of pair-wise potentials yet much shorter than the scattering length. Such systems can display universal 
characteristics independent of the details of the short range interactions. In particular, the energy per particle 
is a fraction £ of the Fermi energy of the free Fermion system. Our main result is that for space dimensions D 
smaller than two and larger than four a specific extension of this problem readily yields £ = 1 for all D < 2 
whereas £ is rigorously non-positive (and potentially vanishing) for all D > 4. We discuss the D — 3 case. A 
particular unjustified recipe suggests £ = 1/2 in D = 3. 

PACS numbers: 



I. INTRODUCTION 



Bose-Einstein condensates (BEC) in dilute atomic gases 
were achieved in 1995 for rubidium, sodium, and lithium [1]. 
Recently, with the production of cold atomic gases close to 
the Feshbach resonance [2], these systems provided a vehicle 
for the study of the BCS to BEC crossover. The (BCS) super- 
conducting and Fermi superfiuid are mixed systems harboring 
quasi-free and pair correlated fermions. In the dilute strong 
coupling molecular BEC, all fermions are relatively tightly 
bound into pairs forming a unique macroscopic state. Several 
theoretical studies have been carried, e.g. [3]. Many of these 
works report on very interesting universal aspects of the BCS 
to BEC transition. In such dilute systems with an inter-particle 
separation much greater than the range of the pair-wise poten- 
tial yet far shorter than the scattering length, the value of the 
energy density at this crossover is independent of specific de- 
tails. This universal energy density is only dimension depen- 
dent. Direct results on the three dimensional system are hard 
to obtain and numerical works have been extremely valuable. 
As is well known from a multitude of other arenas, a dimen- 
sional generalization of original problems posed in three di- 
mensions to arbitrary dimensions often allows us an analytical 
access to original three dimensional problems (e.g. the well 
known e = 4 — D expansion, with D the spatial dimension, 
which has been extremely fruitful in statistical and quantum 
field theories, e.g. [4]). Here, we follow suite and couch the 
BCS to BEC crossover problem in arbitrary dimensions (D). 
We report new results on the energy per particle at the on- 
set of this crossover in such an arbitrary dimension (suitably 
defined by the appearance of zero-energy two particle bound 
states). This extension enables us to examine crossovers in 
low dimensions where no condensed phases occur due to the 
enhanced role of low energy fluctuations. We will illustrate, 
both by exact variational bounds and normalization consider- 
ations, that in all dimensions D > 4 the energy per particle at 



the onset of the transition is rigorously bounded from above 
by zero. By contrast, due to localization tendencies in low 
dimensions (D < 2), each particle carries, on average, the 
mean energy of a free Fermi system at the appearance of the 
first two particle bound states. The physically pertinent well 
known question [5] concerns the energy per particle in actual 
transitions occurring in D = 3 dimensions. The determina- 
tion of the energy per particle here is far more difficult. In 
D = 3, we depart from the more rigorous results in higher 
and lower dimensions and present a heuristic argument for 
recently reported numerical results [6]. These independent 
heuristic arguments bolster the the result attained by dimen- 
sional interpolation. The average of the exact bounds on the 
energies in D = 2 and D = 4, which is half the free fermion 
energy per particle, is not far removed from the numerical re- 
sults reported in three dimensions [6]. In [6] the fraction of a 
half, derived here by (i) dimensional interpolation and (ii) an 
independent heuristic argument, is replaced by £ w 0.44. 
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II. OUTLINE 



We begin, in section(III), by formulating the problem and 
briefly reviewing numerical results. Next, in section(IV), we 
turn to an exact variational bound concerning this transition in 
high dimensions and an easier result relating to bound states 
in low dimensions. Finally, in section(V), we present heuris- 
tic arguments for D = 3. These further correlate the average 
energy per particle with the observed differences seen in finite 
size system with an even or odd particle number. In a brief 
appendix, we provide the specifics of the exact solution of the 
zero energy bound state problem in a D-dimensional spheri- 
cal potential well. Throughout the text, in order to make the 
physics and scaling very transparent, we will often use sim- 
pler forms. Nevertheless, at the end of all calculations, we 
will demonstrate that our results go unchanged with the inser- 
tion of the exact zero-energy bound state solution. 



III. FORMULATING THE PROBLEM: INTRODUCTORY 
AND GENERAL COMMENTS 



We start by writing down the general Hamiltonian. describ- 
ing N — 2n spin- 1/2 Fermions in a box of size V = L 3 , 
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Here, /, J (or i, j) run from 1 to (or n) respectively. The 
phase space coordinates r j , Pi and , are the positions 
and momenta of the "spin up" and "spin down" atoms respec- 
tively and \—gV{vi : Y'j)] is the attractive short range potential 
in mixed pairs. These conventions for the coordinates [using 
upper case characters to describe the total system (both spin 
up and spin down) quantities and the use of lower case charac- 
ters (either "primed" or "un-primed") for spin-down and spin- 
up coordinates] will be consistently employed throughout this 
work. 

Atoms of identical spin polarization cannot be in a relative 
S wave state where the short range interactions operate and 
thus, in Eq.(l), we set V(rj,rj) = V(r^.,rJ) = 0. The cou- 
pling g is tuned to g* where the first (zero energy, S wave) 
two-body bound state appears and the scattering length a di- 
verges. An illustration of the standard two-particle zero en- 
ergy bound state in three dimensions is provided in Fig.(l). 
With this identification of g* at hand [5], the BCS to EEC 
crossover problem is now formally defined in arbitrary dimen- 
sion D, although, due to the enhanced role of low energy fluc- 
tuations, actual condensates may form only in D > 2. Our 
physical interest is in D = 3. To make g dimensionless in 
Eq.(l, we scale V by (mrg) -1 with r$ the range of the poten- 
tial and to the particle mass. We assume a dilute system with 
inter-particle separation d = 1/kp, much larger than ro. In 
practical terms, this implies that we will always consider the 
limit tq — > 0. To avoid the well known formal "collapse prob- 
lem" in which all particles may sit within the attractive poten- 
tial well [as the attractive V{ri,r' j) < for |r,— r'j\ < ro, we 
may envision a state in which all particles sit within a sphere 
of radius ro], leading to a divergent negative energy density, 
the problem is formulated with the provision that the limit 
ro — > is taken before the thermodynamic limit is considered 
[5, 7]. With these definitions in tow, the problem emulates 
what transpires in many dilute fermionic systems (cold atom 
traps, neutron star layers) with separations far larger than the 
range of the pair potentials yet far smaller than the scattering 
length. 

In three dimensions, scaling arguments then imply a uni- 
versal form for the energy per particle 

E 
N 

with £p the Fermi energy of the free fermi system. [8] The 
fraction £ is independent of the specific potential V(r) chosen 
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FIG. 1: A schematic (thick lines) of the binding spherically symmet- 
ric potential V^r, — r'j) = V^(ry) (in three dimensions) which 
is non-zero only inside a sphere of radius ro about the origin, 
V(r) = — Vo@(ro — |r|). A cartoon of a scaled very weakly bound 
s-wave state (x( r ij) — \ r ij \4>(. r ij )) is shown. To attain exactly one 
zero energy bound state we need to fit exactly a quarter period for 
x(|r|) ~ sin fir and set uro — ^ with k = \/2pVo/h, Here, 
fi — m/2 is the reduced mass of a fermionic pair. The zero-energy 
bound state wavefunction in general dimension D is derived in the 
appendix. 



and is the same for all short range potentials which have a zero 
energy S-wave bound state. Cold dilute Fermionic atoms dis- 
play many fascinating universal (and not universal) features; 
see, e.g., [3]. Here we will focus just on the above £ parame- 
ter. 

In actual traps, an excited zero energy state V'cK 1 ') P la Y s a 
key role. The continuum states are effected by the extended 
zero energy state and have negligible overlap with lower, 
tightly bound, states of size ~ r which can be ignored un- 
less we consider very long time scales. As any zero energy S 
wave bound state has the same form ipo(r) = ipo(r) = A/\r\ 
outside the range of the potential for |r| > ro, we assume just 
one zero energy state. For g = 0, there are two free species 
with n = N/2 fermions in each, and relative to Eq.(2), with 
the standard results 
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where pp and Sf = Spr are the free particle Fermi momentum 
and Fermi energy respectively. Throughout, we setfcl. Thus, 
by definition, 



£ = 1 for g = 0. 



(4) 



For small g > (weak attractive potentials), the BCS 
wave function is adequate. By contrast, for strong coupling 
<7|V(|r| < ro ) | >> Sf, up-down spin pairs tightly bind into 
dimers of sizes of order, ro, the range of the potential V. The 
n = N/2 dimers in the deepest bound state then behave as 
point-like bosons and undergo BEC (Bose-Einstein Conden- 
sation) to the traps' ground state with Pdimer ~ so that 
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with B.E. the binding energy. 

The focus of the current work is on the crossover between 
the BCS and BEC regimes wherein the initial Cooper pairs 
become tighter dimers [9]. The existence of a bound state ren- 
ders a weak-coupling perturbation series inappropriate. Sim- 
ilarly, strong coupling/tight binding schemes also fail. As g 
decreases so does the binding energy. Once \B.E.\ < £p, the 
size of the bound state becomes larger than the average inter- 
particle distance: 0(7nep)~ 1 ^ 2 ~ d and antisymmetrization 
of spin up (and of spin down) atoms — which is negligible for 
small dimers — raises E/N above = 0. 

The difficulty and interest of the problem motivated sev- 
eral calculations. A semi-analytic approach utilizing the 
Pade approximation has been discussed at length in Ref. 
[7]. Numerical Fixed Nodal Planes-Green Function-Monte 
Carlo (FNPGFMC or MC) calculations were used to estimate 
E(N)[6]. The ground state wave function of N fermions in 
a periodic box of size L was estimated by evolving (in imag- 
inary time t = —it) an initial trial function ^(rj, r'j; r = 0) 
with cx-p(—Hr) . The short range potential Vfa — r'j) was 
chosen to have precisely one zero energy bound state. The 
lowest energy obtained to date via these methods corresponds 
to £ = 0.44 in Eq.(2) which is consistent with present im- 
precise experimental values. The L/d ~ n 1 / 3 = 80 1 / 3 ~ 3.5 
used there implies a box size L smaller than the (infinite) scat- 
tering length. This still allows for an extraction of reliable 
data. 

In general, the ground state energy E/N decreases contin- 
uously with increasing g. The instability against the forma- 
tion of Cooper pairs for any attractive potential adds a singu- 
lar ~ exp(— g c /g) (with g c a constant) term to E/N with the 
continuity in g maintained. 



IV. EXTENSION TO ARBITRARY SPATIAL 
DIMENSIONALITY D 

While the original problem is posed in D = 3 dimensions, 
it is instructive to address it in any (continuous) D with the 
simple extension of Eqs.(2,3), 

E D 

n = ^WTY) £F - (6) 

To set the ground for future notation, we mention that the 
simple D dimensional extension of the latter part of Eq.(3) 
has the total fermion number (per spin flavor) within a Fermi 
sphere of radius pp given by n — Vtp,L D 'pp / '(2ir) D . Here, 
fljj = 7: d / 2 /T(y + 1) is the volume of the D dimensional 
unit sphere. 

In what follows, we (i) establish by localization tendencies 
in dimensions D < 2 that < 2) = £(g = 0) = 1. We 
then (ii) illustrate that £(£) > 4) < by relying on diver- 
gence, at short length scales, of the wave function normaliza- 
tion. This simple scaling consideration is then made rigorous 
for all continuous dimensions D > 4 by employing a varia- 
tional bound with the use of an "orbital" Slater determinant 
wavefunction. (Regretfully, we have not been able to devise 



a wavefunction which will lead to new stringent variational 
bounds on the three-dimensional problem.) 

Putting all of the pieces together, we find by extending the 
problem to any dimension D while insisting on having pre- 
cisely one zero energy bound state that £ is non-trivial only 
within the interval 2 < D < 4. From this perspective, the 
only interesting integer dimensionality is, in fact, the original 
D=3! 

There may be other extensions of the problem which might 
be non-trivial in all dimensions. Many systems simplify at 
D — > co and a mean field approach may apply there also to 
Fermionic problems [12]. It has been argued that a certain 
class of diagrams dominates in this limit the perturbative se- 
ries in an appropriate effective field theory and can be analyt- 
ically summed in certain fashions to yield £(£) = oo) = | 
[13] or | [14]. 

In what briefly follows, we employ the above continuation 
with a primary focus on energetics (i.e., the existence of a 
zero energy S wave bound state) rather than on the infinite 
scattering length. We find the resulting simple values £ = 1 
for all (continuous) dimensions smaller than D = 2 and £ < 
for all D > 4 interesting on their own right. 

A. Low Dimensions (D < 2) 

Localization tendencies present in low dimensions afford 
us with a direct result. Purely attractive potentials V (r) have 
zero energy bound states in one and two dimensions [10], for 
any strength and g. It follows that a bound state appears imme- 
diately as the two body attraction is introduced, i.e. g* = + . 
As the energy fraction £(g) is continuous in any dimension D, 
we find that 

e = £(.9 = 0) = l, D<2. (7) 

Similar incarnations of localization tendencies in low dimen- 
sions have proven very useful in other arenas (e.g. [15] in 
which a scaling theory of localization suggested D = 2 as the 
marginal dimension for the appearance of metallic states). 

B. High Dimensions (D > 4) 

1. Energy attained by Normalization conditions 

Next, let us assume that the potential V vanishes outside 
the range ro- Volume normalization factors aside, the spin-up 
spin-down, zero energy, pair wave function for relative sep- 
aration |r| > 7'o is 4>(r) = njrif=2) ■ with Ay a dimension 
dependent numerical constant. The following is a very inter- 
esting observation: the integral representing the normalization 
of this zero energy wave function, [J d D r\<f)(r)\ 2 ] diverges at 
small r for D > 4. We consequently find that in D > 4 di- 
mensions most of the normalization of the zero energy bound 
state wave function is concentrated near |r| = r$ and the 
dimers formed are actually compact — just like for large cou- 
pling in D=3. Hence, we can again neglect the overlap of 
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the different dimers and attendant kinetic energy due to anti- 
symmetrization in the spin up coordinates {r^} and separately 
in the spin down positions {r^}. The strong coupling result 

Eq. (5) applies and we have: E/N — — < 0. Thus, we 
find that £ may vanish (or become negative) for all dimensions 
greater than four, 



£ < , D > 4. 



(8) 



In the marginal dimension D = 4, the size of the dimers is 
just barely ro (where the attractive potentials operate) and the 
upper bound may be more easily saturated with a potentially 
vanishing £(£) = 4). 



2. A Variational Upper Bound 

To substantiate the above arguments, and directly derive 
Eq.(8), we next employ a 2n-body ansatz wave function 
^(ri, r 2 , r,, . . . r„; r[, . . . , y' 3 . . . r' n ) and show that in all 
continuous dimensions D > 4, we have H\^f) = 0. The 
ansatz has the form of a Slater "orbital" determinant. It rep- 
resents the state of n dimers, say, <p(ri — r^) . . . <j>(r n — r' n ) 
with 4>{\Ti — r'j |) the zero energy bound state, antisymmetrized 
over all the n! permutations: P = i — > (p(i)), i — 1, ... n of 
the n spin-up atoms: 



^(-1) P $(F); *(P) = [] ^i-t' p{i) ) (8) 



with (— l) p the parity of the permutation P. Here, ^(r^ — 
r p(i)) represent normalized, zero energy, pair wavefunctions. 

The further required antisymmetrization over the n spin- 
down atoms simply scales ^: ^ — > [(n!)^] and is redundant. 
We note that ^ can be expressed as a determinant, 
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We can use the variational ansatz in any dimension D to 
obtain a variational upper bound on the N=2n-body ground 
state energy and on the corresponding £: 



Clearly, 
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(11) 



with a similar expression for (^|^). 

The local potential limit, say, r - 
spherical potential 



in the "square well" 



V(r) = -Vb6(r - |i 



implies for all dimensions D > 2 the following "selection 
rule" for non-vanishing ($>(P')\Vij\Q(P)) matrix elements 
appearing when the 2n-body Hamiltonian (1) is substituted in 
Eq. (11) above. 

Before going into the details of the variational calculations, 
let us quickly give the reader a glimpse of the final results 
derived below and their implication. The calculation detailed 
in this subsection illustrates that for the value of Vq in Eq.(12) 
(the analog of g in Eq.(l)) which secures a single zero-energy 
bound, state we have a variational state (i.e. of Eq.(9)) for 
which the corresponding energy E var = 0. This result then 
implies (by Eqs.(6,10)) that £ < 0. 

For any Vi.j among the n 2 spin-up-spin-down potentials 
we have, for all D > 2 with a local V(r) producing one zero 
energy bound state, the following relation 

($(P')\Vi j|$(P)) = if nop(i) = 7 and no y'(i) = j. 

(13) 

To prove this claim, we perform first (among the 2n 
dP Y^d D y[ integrations involved in evaluating the above ma- 
trix element) the integrals J dP 1 r< J dP r'j over the coordi- 
nates Yi and y'j appearing in the particular potential term 
V(yi — r^) considered. 

Both Yi and y'j appear in the product of Eq.(8) but, by 
assumption, in different <fi factors. The same holds for $(/"). 
Thus, the integration over and r'j is of the form, 



rf 3 r;[0(r, 



p'(i) 
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/)], d4) 



with p^ 1 the permutation inverse to p. The condition for hav- 
ing one zero energy S wave bound state in the square well 
potential of (12) in the D=3 case is well known (see Fig.(l)), 



Vq 



(15) 



By dimensional considerations, the same condition for one 
zero energy bound state holds in all dimensions D > 2 
wherein the constant 7r 2 /4 is replaced by some other, dimen- 
sion dependent, numerical factor cd, 
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CD 



(16) 



(12) 



In the general D dimensional problem, cd of Eq.(16) is fixed 
by normalization of the zero-energy wavefunction </>(r). In the 
appendix, we briefly present the solution to the a zero-energy 
bound state in a D-dimensional spherical potential well. 

By changing variables (r,,rj-) — > (r^r^) with r^j = 
(Yi — Yj), the local Vfaj) implies that r$ = y'j in the argu- 
ments of all the four <fi factors appearing above which gener- 
ically keep finite arguments and bound cf> values. Neglect- 
ing small variations of </>(r) away from r ~ 0, the integral 
J d 3 Y V(y) over the D dimensional sphere of radius ro where 
the square well potential is non-vanishing yields the factor 
VqtEVId- This multiplicative factor is, by virtue of Eq.(16), 
proportional to [crjr Q D ]. In dimensions D > 2, this fac- 
tor vanishes as ro — > and, as claimed earlier, so does the 
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complete matrix element. This is not the case in dimensions 
D < 2 for which £ = 1. 

Next, we consider dimensions D > 4. In this case, the zero 
energy bound states 

tin, r' pW ) = iT D/2 </)(r ? - r' p(l) ), (17) 

with (see the Appendix for details), 

d>(r i -r> pii) ) = A > \r i -r' p{{) \-( D -V 
for |r| = |r ; - r' p(i) | > r ; 
and ^( ri -r; w )~A>(r - (i3 - 2) ) 

forO<|r|<r , (18) 

where A> ps ( r o )■ Exact forms are provided by 

Eqs.(A3,A4, A8, A9). The comparison between Eq.(18) and 
the exact forms is provided in Eq.(A7). The scaling of the 
results with ro (which we will shortly obtain) becomes more 
transparent with the use of Eq.(18). With the incorporation of 
Eq.(A7) and its ensuing discussion, the results which we will 
obtain using Eq.(18) will further enable as a rigorous upper 
bound on the matrix value elements to be discussed. As seen 
from Eq.(18), the wavefunctions <fi are strongly localized at 
r ~ ro implying stronger selection rules for non-zero matrix 
elements: 

($(P')|$(P)) = unless P = P', (19) 

and 

($(P')\V(ri - r^)|$(P)) = unless j = p(i)=p'(i). 

(20) 

Eq.(20) implies that also (<I'|iJ|^I'}, the numerator of 
Eq.(10), has only diagonal contributions. Namely, 

<$(P')|P|$(P)) = unless P = P' . (21) 

To illustrate Eq.(19), let us assume that P and P' differ 
minimally: p(i) = p'(i) for all i > 2, butp(l) = 1, p(2) = 2 
andp'(l) = 2, p'(2) = 1. The overlap in Eq. (19) is then 

L- 2D J d D v 1 d D v' l d D v 2 d D v' 2 4 ) {v 1 - r;)0(r 2 - r' 2 ) 

x^(n -y' 2 )4>(y 2 - t' x ) 

L- nD [] / d D v t d D v' 3 m^))m,p{i)))- (22) 

I 



where we invoke Eq.(A4). In particular, as seen from the ap- 

1/2 

pendix, the wavenumber k = fro- In Eq.(23), we employ, 



The last line of Eq.(22) denotes a normalized integral (equal 
to one). In the first integral over the 4 variables {ri 2, r'1,2}, 
the number of integration variables is equal to the number of 
wavefunctions {<^(r; — r'j)} appearing in the integrand. For 
each of the two particle wavefunctions <^>(rj — r'j), we in- 
sert Eq.(18) (or the exact Eqs.(A3,A4) derived within the ap- 
pendix) with Ay ps Tq 1 . Insofar as scaling with ro is 
concerned, all integrals over products of wavefunction forms 
of Eq.(18) multiplying A> amount to product of either in- 
tegrals of the type /> = Ji r i >ro d D r\r\ 2 ~ D or of the form 

I < = rir D f, , . dPv. Constant factors aside, both 7> and 

^ U J|r|<ro ^ 

I < scale as r 2 .. The four normalization prefactors of A > (with 

(Y) /2~) 2 

Ay ps Tq ' ; ), originating from the four factors of cf> in 
the top two lines of Eq.(22), lead to an additional factor of 
rj) 4 ' in tow. Thus, in the final analysis, when the integrals 
are segregated into all possible terms for |rj — r'j\ > tq and 
for |rj — r'j j < ro and Eq.(18) or Eqs.(A3, A4)) are inserted, 
we find that all terms scale as A ps rg D — > as ro —> 0. This 
signals that the factor in the first two lines of Eq.(22) vanishes 
as ?'o — > and thus so does the entire overlap of Eq.(19). We 
can readily verify that for larger mismatches between P and 
P' the overlap vanishes as a higher power of ro. Now, in the 
diagonal case, with p'(l) = 1 and p'(2) = 2 the product of 
the four normalization factor cancels and ($(P)|$(P)) = 1. 
[16] 



Turning to Eq.(20), we assume that p(i) = j but p'(i) 7^ j 
and prove that the matrix element (P'|Vij|P) vanishes for 
D > 4. Recall that when also p(i) ^ j we showed that this 
matrix element vanishes in all D > 2. Unlike that previous 
case, we have here only three (rather than four) <fi functions in 
which r ; and/or r'- appear. Factors of L aside, the relevant two 
integrations on r; and rj in the analog of Eq.(14) are now: 



(23) 



once again, p 1 for the inverse permutation, r = r; — r'j is the 
argument of the square well potential, and 0(r; — r'j) = 4>{r). 



d D r, / d D r' 3 Vir, - r^fa - r' p{l) )<f>(r p ,- Hj) - r'^r, - r'j) 



A^o I d D n j d D rQ(ro - M^Kr) 1 "-* J f ^M]^ - r^)^-,^ - r^), 

I 
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In what briefly follows, to flesh out more crisply the simple 
scaling form of this term, we may replace the exact Bessel 
function form of Eq.(A4) and instead invoke Eq.(18) to re- 
place the last factor (4>(r)) for r < tq in the upper line of 
Eq.(23). The result which we will obtain in this fashion will 
go unchanged if we employ the exact Bessel function form 
of Eq.(23), see [17]. With the insertion of Eq.(18), the d D r 
integration of Eq.(23) yields 

V r d D r [4- 2 r { t D) ] = n D V Q rlrf- 2 



ft2 o — 

— U d cd r 

171 



(for D > 4) , (24) 



in the limit ro — > 0. Here, the definition of cd was invoked 
from Eq.(16). The peculiar form given in the square braces 
of Eq.(24) follows from Eq.(18) with the approximate A > m 
(t-q ). The incorporation of the exact Bessel function 

form of Eq.(A4) in Eq.(23) whose derivation is detailed in the 
appendix does not alter this simple result [17]. This concludes 
the proof of Eq.(20). 



With all shown thus far, vanishing in Eq.(23) may be 
avoided only when p'(i) = j = p(i). In this case, the 
and r'j integrations involve only two identical cf> functions and 
become (invoking Eq.(A4) once again), 

/^r^r0 2 (r)Vb-e(ro-|r|) 
= V A% [ d D r t [ d D r [{Krf-^Jo^r)] 2 

J J0<\r\<r a 2 

Vr;Fo-(25) 



By virtue of the zero-energy Schrodinger equation that 
</>(r) satisfy, the contribution of Eq.(25) (albeit divergent) 
identically cancels against the kinetic terms Ti + T- = 
— (ft 2 /2m)(V 2 + Vj 2 ). The coordinates Vi and r^- now ap- 
pear in both Q(P) and $'(P) only in the above two identical 
<f> functions on the left and on the right, Ti + T'- can be paired 
with the Vij considered here into Hij = Tj + Tj + Vij. In 
the above shorthand notation, the matrix elements of Hij are 



J 



mP')\H^(P)) = (<t>{l,p'{l))...<p{n,p'{n)) | V i% 



r 



2m 



(V 2 



Vf) U(l,p(l))...0(n,p(n))). 



(26) 



The Schrodinger equation for 4> with the assumed small bind- 
ing energy (eventually e — > 0) which defines our problem is: 



This holds for each $(P) and, hence, for the full Slater orbital 
ansatz of Eq.(9) ^ in D > 4 dimensions, 



2m 



+ V(r i -r' i )}<P(\r i -r' j \) = -e<j>(\r i -r' j \), 

(27) 

or Hijtpij = —e<f>i.j. Consequently, the matrix element of 
Hij between the two <pij in and in $(P') is just a num- 
ber (e). Omitting, then, <pij = <fi(i,p(i)) = <fi(i,p'(i)) from 
the product of n functions <j> in $(P) and $(P') leaves an 
overlap of two monomials each composed of a product of the 
remaining (n — 1) orbital wavefunctions <f> in both the bra and 
in the ket, i.e. ($(P^_ 1 )|$(P„_i)} with P„_i and P' n _ x , the 
permutations of the remaining [n — 1) elements obtained by 
removing i — » p(i) = j and i — > p'(i) = j. Eq.(19) implies 
that this overlap vanishes unless these remaining permutations 
are identical. Inserting i — > p(i) = j and i — > p'(i) = j im- 
plies that P = P'. It follows that in the numerator of Eq.(10) 
we have, similar to the denominator, only diagonal terms. 

By definition, the permutations (i,p(i)) reproduce the n 
distinct (i, j) terms for given i. We may, therefore, operate on 
$(P)) on the right with the n potential terms V(i,j) and the 
?i kinetic energy operators T^j with matching <f>(i,p(i)) terms 
in $(P). In the process we employ each of the Laplacian 
terms ({V 2 } and {V^ 2 }) exactly once. All further unmatched 
V(i,j) have vanishing matrix elements. 

Hence, 



H\V) = E(2n) 



-ne|*) 



(29) 



Thus, the total energy of the system is just that of n inde- 
pendent Bose-Einstein condensed dimers. When the binding 
energy e — > we have E{2n) — > and for the variational or- 
bital Slater determinant of Eq.(9), we find £, var {D > 4) = 
as claimed. By the standard variational theorem, this proves 
that that for continuous D, 



£{D > 4) < 0. 



(30) 



P|$(P)) = -ne|$(P)). 



(28) 



The exact variational result fortifies the conclusion arrived at 
earlier (subsection(IV B 1), with Eq.(8) in particular) via nor- 
malization considerations. In the final analysis, normalization 
was present in our slightly elaborate variational bound and ul- 
timately dictated the scaling form of our expressions with vq 
(see, e.g., Eq.(22), [16]) as ro — ► 0. 



V. HEURISTIC CONSIDERATIONS FOR £(D = 3) 

The extension of the problem to arbitrary continuous di- 
mension D allowed us to firmly establish that £ = 1 in D < 2 
and £ = in D > 4 (with the particular point D = 4 accessed 
by general normalization conditions; our rigorous variational 
bound £(£) > 4) < held for all continuous D > 4). These 
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results suggest that in a lowest order e expansion interpolating 
between these two limits, £(£) = 3) « 1/2. [Or, if the up- 
per bound in the variational inequality for all continuous large 
D, £(D > 4) < is not saturated for the marginal D = 4, 
that, more generally, £(D = 3) < 1/2.] This is, remarkably, 
not far from the results obtained by very intensive numerical 
work [6] (which led to an upper bound of £(£) = 3) « 0.44). 
Unfortunately, we have not been able to attain direct potent 
bounds on the three dimensional problem which are close to 
the numerically attained values. In the current section, we at- 
tack the harder D = 3 problem by a heuristic approach. This 
approach, although imprecise, may shed light on some of the 
observed physics in three dimensions. Similar to the simple 
minded dimensional interpolation result it, too, suggests that 
£(D = 3) w 1/2. As an additional bonus, this allows for a 
heuristic argument for even-odd oscillations numerically ob- 
served. 

Unlike the normalization of the zero energy bound state, 
which for |r| > ro is uniform in |r| = |r» — r'-|, the kinetic 
energy, 

J d 3 r(V0(r)) 2 = J d 3 rj^, (31) 

is mostly concentrated at small |r| ~ ro << d values. In the 
last equality of Eq.(31), we inserted the well known (D = 3) 
zero-energy S-wave wavefunction form outside a potential 
well, <fi = 4^r. As a consequence of this concentration of ki- 
netic energy, a zero energy bound state of a mixed (i,f) pair, 
can manifest over any smooth background function B of all 
other variables by a nrz~p enhancement at distances — r'j | 

smaller than typical particle spacing. (Due to antisymmetry, 
the i th spin up atom is equally likely to bind with any one of 
the n down spin atoms at r^, with j = 1, n.) 

This is indeed the case for the wave functions serving as 
the staring point of the above MC where correlations in mixed 
pairs are introduced via a Jastrow product J, 

^=BxJ = B H tdn-r'jl), (32) 
3=1, ■■n 

where the background B of free Fermion Slater determinants 
or BCS form accounts for the antisymmetrization. 

The other (n — 1) factors of the form , — 2_t with k differ- 

ent from j still allow for "J/ ~ ^rz~p | f° r I r i — r j I < < d. Thus, 
let the remaining (n — 1) down spin atoms form a cubic lattice 
of spacing d, the triply periodic product, Ilj = i n -\ <Ml r j — 
rj|), is then invariant under the point symmetry group of the 
simple cubic lattice. Its finite extrema are at the centers of the 
n cubes (each of volume d 3 ) and C varies slowly for |r| << d. 

This suggests (yet does not justify) the following ad-hoc 
approach which yields £ = 1/2. The potentials V(|rj — r^|) 
involve atom pairs yet the kinetic energy is a sum over single 
atoms: 

p? = E p? + E (pP! 

^ 2m ^ 2m ^ 2m 

Z=l,..JV=2n i=l, ...n j=l,...n 



We formally associate kinetic energies with pairs by using the 
identity : 

Ei^S^TT) £ (p '- p -' )2 ' <34) 

\ > I<J=1,2...N 

where P| o( = (J^Pi) 2 = was subtracted cancelling all 
I < J terms (p/ • pj). Next, we separate the contribution to 
the Hamiltonian of same (up-up, down-down) and mixed spin 
pairs with the potentials included in the second term: 

/ } T ~ -i ; J2 (Pi - Pi) 2 

v ' %<l=\,...n 

+ 2m(N + l) E (Pi-Pi)' 

v ' j<k=l,...n 

For g = the kinetic terms in H same and H mixe( i are equal 
up to 1/N corrections. 

Next, we note that if the following assumptions are made: 

(i) Just as for isolated up down pairs also in the N- 
body ground state "J, the attractions between different 
atoms cancel the kinetic energy of the relative motion, i.e., 
{^\H mzxed \^) = 0, and that 

(ii) turning on the coupling g does not change the expecta- 
tion value of Hsame, then: 

{E/N)\ N ^ en = ^l\e F 

+ 06) 

For odd N, say n+ 1 spin up and n spin down atoms, only n = 
N/2 — 1/2 bound pairs can form The decrease of potential 
attraction energy by (N — 1) /N produces a "gap" 

E E 

-^-|AT=2n+l ~ -j-j\N=2n ~ £F, (37) 

consistent with the finding of the MC calculations. Further- 
more, 

£ 1/2 (as N oo), (38) 

which is barely consistent with the variational bound £ < 
0.43 — 0.45 from the above calculations. [18] 



VI. CONCLUSIONS 

In conclusion, we examined an extension of the BCS to 
BEC crossover problem to arbitrary dimension D. We re- 
port new results on the ground state energy per particle at 
the crossover point by employing direct variational bounds in 
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high dimensions and by considering the consequences of lo- 
calization in low dimensions. In particular, we find that the 
ground state energy per particle at the onset of the BCS to 
BEC crossover is zero (or negative) in all dimensions D > 4 
while it is the energy of a free Fermi system in all dimensions 
D < 2. The interpolation of these bounds to the physical three 
dimensional problem leads to an energy per particle which is 
half that of the free fermion energy and is close to current nu- 
merical results. We outlined a simple heuristic argument for 
the same result which further mandates even-odd variations 
which are indeed seen numerically. 
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with n and the constant prefactors A > and A < determined by 
continuity at r = r and global normalization. Here, Jd__ 1 

is a Bessel function of order (y — 1). (Inserting Ji(nr) = 

[(ftr) -1 / 2 sin Kr], the pertinent three dimensional result, illus- 
trated in Fig.(l), is recovered.) Continuity at r = tq restrains 



(r < r ) « A 



-2--D 



>'0 



(A5) 



For small r (i.e. r — > 0), employing the asymptotic form of 
the Bessel functions, we find that 



(A6) 



Anywhere within the potential well (r < r ), the wave func- 
tion 



APPENDIX A: THE ZERO ENERGY BOUND STATE 
PROBLEM IN A D-DIMENSIONAL SPHERICAL 
POTENTIAL WELL 

Here, we provide the solution of the zero energy bound state 
problem in D dimensions for the spherical potential well of 
Eq.(12) wherein the depth of the well (Vo) is tuned to get a 
zero-energy bound state. The solution is a straightforward 
extension of the standard D = 1,2,3 dimensional spherical 
potential well problems. In these, without any angular de- 
pendence, in the the "S-wave" representation (i.e. the scalar 
("I = 0") representation of the SO(D) rotation group), the 
D-dimensional Laplacian is V 2 = {-^ + ^r-jj^]- This ap- 
pendix explicitly illustrates how the numerical constant cd of 
Eq.(16) may be determined (from the implicit Eq.(A8)). Its 
results further allow us to compare the approximate form of 
Eq.(18) introduced within the text to allow a clear understand- 
ing of the scaling form of the overlap integrals with the exact 
form of <p(r) (Eqs.(A3,A4)). In the aftermath, it will be shown 
that the incorporation of the exact functional from does not 
change the scaling results derived in the main text (following 
Eq.(18)). 

We proceed with the solution of the translationally invariant 
problem specified by the spherical potential well of Eq.(12). 
By translational invariance of the center of mass, the two body 
wavefunction 



(x,x') = L- D/2 < 



(Al) 



with r = x' x. In the potential-free region (r = |r| > 
ro), the wavefunction (j> satisfies a single particle Schrodinger 
equation with a reduced mass [i = m/2, 



d 2 ,. . D-l d ,, , 



K^{r) = 0, (A2) 



with k = 0. Within the spherical potential well (r < ro), we 
have Eq.(A2) with k 2 = (mVo)/h 2 . We immediately find that 



(r > r ) = A y 



„2-D 



while within the potential well 

4>{r < r ) = A<[( K r) 1 "* Jd^kt)], 



(A3) 
(A4) 



r 

'ro' 



(r < r ) = f(-)A >r 2 a - D , 



(A7) 



with /(0 < z < 1) a bounded function satisfying /(l) = 1 



and/(0) 



r D-2 
lfl_ 



Within the text following Eq.(18), 



l > r{§ )2i 

we set / to be one and we attained vanishing overlaps in 
Eq.(24) in the limit ro — » 0. The incorporation of the more 
explicit, yet bounded, f(z) given by Eq.(A4) does not alter 
this result. Formally, our results derived in the text can be 
made rigorous by deriving an upper bound which differs from 
the results given in Eq.(24) by a factor of |/| ma x, the maximal 
value of |/| . As Eq.(24) vanishes by simple scaling, the incor- 
poration of a finite factor leads to a vanishing upper bound, 
proving that Eq.(23) vanishes in the limit of small ro. Al- 
though inconsequential for this bound (all that matters is its 
bounded norm), the prefactor /(r/ro) is everywhere positive 
as the zero energy ground state is nodeless and of uniform 
sign. 

For completeness, we now outline the general solution. 
Matching the two functional forms of <j>(r) and their deriva- 
tives at r = ro (and simplifying with the aid of the standard 
Bessel function recursion relations) leads to the implicit equa- 
tion 



(1 - y)J^_!(Kr ) 
:(nr )[JD («x ) + Jd_ 2 (kxo)]. 



(A8) 



The solution of Eq.(A8) enables a determination of k and 
thus of the spherical potential depth Vb and the constant co 
in Eq.(16) which ensures a zero-energy "S-wave" state in D 
dimensions. The amplitudes A > and A< are then determined 
by 



„2-Z5 



(kxq) 1 ^ JD_ 1 (nr ) 



(A9) 



in conjunction with the global normalization of (f>(r). Normal- 
ization demands that, up to dimension dependent numerical 



constants, A ; 



D/2-2 



and consequently A < 



-D/2 
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